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Abstract 

The effect of nonlinear self-interaction can be associated with su- 
perluminal velocity of neutrino. The power energy spectrum E = 
p + Cp a is derived from the nonlinear Dirac equation when interac- 
tion term V = X^^f^tp^tj:)' 1 is added to the Lagrangian of a free 
spin-1/2 particle. The superluminal velocity recorded by the OPERA 
and MINOS collaborations is achieved when the coupling constants 
are taken in the range a = 0.4 -j- 1.18 and A = — (0.5 -r 1.6) x 10 -4 . 
The self-interaction Lagrangian V = Xifj-f^tp'ip'-f^'ijj with the coupling 
constant A = - (0.7 4- 0.9) x 10" 4 yields the same result. Scalar inter- 
action V = A( , 0V ; ) 6 an d scalar-vector interaction A(^^) b+1 / (V'V') 6 
cannot be responsible for the observed superluminal neutrino. 



1 Introduction 

Neutrino was believed to be a massless spin-1/2 fermion with energy 

E = pc (1) 

and group velocity v = c equal to the speed of light c = 1 (in relativistic 
units). The modern theory expects, however, that neutrino has finite mass 

H! 

m = m„ < 0.28 eV (2) 
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that implies deviation from the energy spectrum ([[]) and velocity 

dE , . . 

Recent experiments of the OPERA Collaboration [2J revealed superluminal 
motion of neutrino with energy E = 17 GeV at the average velocity 

v = 1 + 2.37xl0" 5 (4) 

while the the MINOS Collaboration detected velocity 

v = 1 + 5.1XKT 5 (5) 

for the low energy neutrino with energy spectrum peaked at approximately 
E = 3 GeV. Superluminal neutrino was also observed in supernova explosion 
SN1987a @|. 

This fact is a serious puzzle to the researchers. Superluminal velocity (j3J) 
cannot belong to a free massive particle with energy spectrum E = \J p 2 + 
whose velocity is always subluminal because 



m? 



1 m 2 1 m 2 , , 

The tachyonic energy spectrum E = ^/p 2 — m 2 results in velocity above the 
speed of light 

1 m 2 1 m 2 , . 

that does not exceed 3 x 10~ 22 even at maximum possible neutrino mass 
m (J2j) and E = 17 GeV. Nevertheless, in the frames of the accuracy of 
measurements, the energy spectrum of superluminal neutrino of OPERA [2] 
and MINOS [3] can be fitted to a power law El E] 

E = p + Cp a (8) 

where the coefficients must be taken in the range [8] 



a 



0.40 -=-1.18 C = 4.15 x 10~ 4 4- 1.5 x 10~ 5 (9) 



and C = (2.2 4-3.03) xl0~ 5 if we choose a = 1. Indeed, neutrino is not 
a free particle, but there are several interesting hypotheses to explain its 
superluminal motion [9]. 
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In the present paper we consider massive neutrino whose Lagrangian 



L = ip(i^d^-m)4> + V (ip,ip) (10) 

includes nonlinear self-interaction term V (ifi, iff\ . There is no additional 
interaction with external fields and the medium, while the superluminal ve- 
locity is hidden in the very nature of neutrino. We need to find the energy 
spectrum of nonlinear Dirac equation 

(<y*^-m)^ + ^ = (11) 

and check whether the Lagrangian (flO|) is adjustable to reproduce the super- 
luminal velocity ©-© detected by the OPERA [2] and MINOS [3j collab- 
orations. 



2 Effect of self-interaction 

Consider the Lagrangian ( flOl) with a simple self-interaction term 

V(^^)=^-^H) b+l ^ = ^7° (12) 



The relevant Dirac equation ffTTj) is written 

(i^d^-m + F)ip = (13) 

where 

F = 7°u; (14) 

and 

uj = \(^i)) b (15) 

The Dirac equation (FTBl at u = has well-known solution in the form of 
plane wave 

= ( t° J exp (-iE t) xo = t^^0o = V™ 2 +P 2 (16) 

Substituting stationary wave function 

if> = ip(f) exp (—iEt) (17) 
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in (ITHjl . we have equation 



(—ia ■ V + (3m) ip — (E + u)ip (18) 

where 

a = M = 7° = -7o (19) 
Substituting a plane-wave bispinor 



<P= ( J exp(ip-r) (20) 



in (fTSj) , we obtain a linear system of equations 

& ' PX = (E + u — m) (ft 
a ■ pcj) = (E + uj + m)x 



(21) 



that has solution if and only if 

E = v 7 ™ 2 +p 2 -co (22) 

where p = \p\ . 

We can estimate the energy spectrum ( l22l) in the frames of mean-field 
approximation if we neglect correlations between the field operators in ( [151) 
and apply effective interaction 

uj ~ = An 6 (23) 

where 



n 



<^V> (24) 

is the particle number density. The latter can be adjusted as n = 1/V for 
1 particle in volume V. For a many-particle system the quantity (|24|) is 
determined according to formula 



n= j2^fJ fkdqk (25) 
o 

where f p is the distribution function in g-dimensional momentum space. The 
latter is the Fermi-Dirac distribution function if it is an ideal gas in equilib- 
rium, however, the OPERA [2] and MINOS j3] collaborations considered a 
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neutrino beam with the energy peaked at k = p, and its distribution function 
corresponds to a delta-function 

h = S (l - ~) (26) 

in 1-dimensional momentum space, so that the particle number density ( 125|) 
is estimated so 

n = £ (27) 

7T 

and, according to (1221) and (1231) . the energy spectrum is 

£ = ^Jp 2 + m 2 - An fe = vV + ™ 2 - 4 P 6 (28) 

3 Vector self-interaction 

Consider more general form of vector self-interaction [TU] 

V = ^T[(^^) b+1 ( 29 ) 
corresponding to the Dirac equation 

{ij^ - m + F) i> = (30) 

where 

F = X (^7^7^)^ (31) 
At b = 1 it is no more than the Heisenberg model of self- interaction [11 J 



F = A 7m (^) (32) 

When we consider 1-dimensional neutrino beam, we take into account 
that gamma-matrices in (l+l)-dimensional representation [10] 



satisfy standard commutation relations 



{ln,lu} = 2v»v V^=(q \ (34) 



Substituting stationary plane wave solution 



u 



ip = exp (ip x x - iEt) p = (p x , 0, 0) (35) 



v 

in the Dirac equation (145]) . we have [10] 

ip x v = (E + u — m) u 
—ip x u = (E + u + m)v 

instead of (1211). while 



(36) 



u = \(\u\ 2 + \v\ 2 ) =A(^V) b (37) 



formally coincides with (TT5j) . The linear system f[36|) has solution if and only 
if condition f[2"2"]) is satisfied. Again, applying the mean-field approximation 

u ~ w„ = A <^ t ^) b = An 6 (38) 

we obtain the same formula for the energy spectrum ([28"]) of a 1-dimensional 
neutrino beam. 

The group velocity is immediately calculated 

" X 4p^ (39) 



a/p 2 + m 2 7T 6 

that tends to 



A6 b _! m 2 

in the ultra-relativistic limit E ~ p ^> m. The latter term in (T4"0"]) is negligible 
even at the upper bound of neutrino mass and superluminal velocity flU- 
(]S]) can be explained by the second term if Xb < 0. Indeed, the power energy 
spectrum (jHj) is compatible with ( 140]) if 

6 = a = 0.4 4 1.18 (41) 

A = -(0.5 4- 1.6) x 10~ 4 (42) 
Particularly, the coupling constant may vary in the range 

A = - (0.7 4 0.9) x 10" 4 (43) 

when a = b = 1 that corresponds to the Heisenberg model ( 132]) . Therefore, 
the origin of superluminal velocity of neutrino [2] [3] can be associated with 
vector self-interaction f ]29]) if the coupling constants are properly identified 

on])-©. 
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4 Scalar self-interaction 

Consider the Lagrangian ( TTUj) with scalar self-interaction 



V 



A 



-{^) b+1 (44) 



6 + 

The Dirac equation (jlip has the form 

(i^dp - m + W) ip = (45) 

where 

W = \(^) b (46) 

Substituting stationary wave function ip = <p(r) exp (—iEt) and the ef- 
fective mass 

m 1f = m- W (47) 

in ( 145]) . we obtain equation 

{-id ■ V + /3m,) cp = Ecp (48) 
that has plane-wave solution 

?=(£)exp(ij?.f) X=^^J> (49) 
for a free particle with the energy spectrum 



£ = VV + ™* (50) 

and relevant group velocity 

p + m*-^ 1 ] (51) 



a/p 2 + m 2 V dp 

Solution (T49l also implies 



# = II0H 2 -||X|| 2 = ^V (52) 

where 

^V= I|0|| 2 +I|X|| 2 (53) 



Applying the mean-field approximation, we neglect correlations in 
and use the effective interaction 

W ^W, = Xn b s (54) 

where 

n s = {rpxp) = -£n (55) 

is the scalar density and n is the particle number density (124"]) . For a many- 
particle system it is determined by formula 

oo 

n -jkS' I m hi,k (56) 



The distribution function /& of a neutrino beam with the energy peaked at 
k = p is presented by a delta-function (1261) in 1-dimensional momentum 
space, so that the scalar density ( ]5"5l) is estimated 



_p m t (p) 

7T E (p) (57) 

Substituting (|54l and ( IBTl) in f j47l) . we obtain self-consistent equation for the 
effective mass 



m* (p) = m — An s = m — A 



p m* (p) 
7T ^/(pa + m 2) 



(58) 



Consider the ultra-relativistic limit E ~ p ^> m*. The energy ( IBOl) is 
reduced to 

m 2 m 2 mW 7 * , N 

£=p+^=p + + — *- 59 

y 2p F 2p p 2p K J 

while the effective interaction (1541 is reduced to 



A 



W, ~ -V (60) 



7T 



and velocity f[5"Tj) is estimated 



2 



« ~ 1 - ^§ (61) 

Although the deviation from the speed of light may be visible at large W*, 
the velocity flBTl) is subluminal at any choice of constants A and b. Therefore, 
the scalar self-interaction cannot be responsible for the superluminal neutrino 
velocity. 
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5 Scalar-vector self-interaction 

Now consider a mixed variant of scalar-vector self-interaction in the form 
[HITS] 

V = X-^—^—j— (62) 
that corresponds to the Dirac equation 



(i^d^ - m + W + 7°w) ip = (63) 

(64) 



where 

u = A {b + 1) ±— } - W = -\b^—\ T 



Substituting a stationary plane wave solution 



ip = ( ^ ) exp (ip- r — iEt) (65) 



and the effective mass (HTI) in (1631) . we have a linear system of equations for 
the spinors 

a ■ px = (E + u - m*) 
a ■ p(f)= (E + u + m*) x 



(66) 



that has solution 



X = (67) 

i? + u + m* 



if and only if the particles has the energy spectrum 



E = y/m\ + p 2 - cu (68) 

where the effective mass is 

= m — W (69) 
According to (I65I) and (1671) we have 

_ II0II 2 + llxll 2 _ V m2 * +p 2 



W> \W-\\x\ 



(70) 
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and the effective mass (|6T?|) is determined by self-consistent equation 



m* = m — Xb 



yj (ml + p 2 ) 



2 ,b+i 



m1 +1 



that is reduced to 



m* = m — Xb 



P 



b+l 



ml +1 



in the ultra-relativistic limit p 3> m*, so that 
~b + 2 



b + l 



-m* — m 



m' = (m, — m) — - 
P 



The energy spectrum (1681) is, then, reduced to 

E^p+^-Xib+l)^ 
2p ml 

At large coupling constant A, when 

m* > m 

equation (1751) implies 



Him* 



6 + 2 p 

and, according to equation (171]) . the neutrino velocity is subluminal 

,2 



t> ~ 1 



m; 
2p 2 



When 



m* — >■ m 



(71) 



(72) 



(73) 



(74) 

(75) 
(76) 

(77) 
(78) 



equation (1751) implies — > and, according to equation (1711) . the neutrino 
velocity is also subluminal 



I! ~ 1 



m 



When 



2p 2 



(79) 
(80) 
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equation (1751) implies 



< ~ — (81) 
V 



and, according to equation (|71|) . the neutrino velocity is superluminal 



2 2 

mf m 



' l! -^2?«2? < 82) 

but it negligible with respect to estimation (J7|) for a free massive neutrino 
and, hence, cannot achieve velocity fll])-© that was observed in experiments 



6 Conclusion 

The model of neutrino with self-interaction can explain superluminal velocity 
dU)-© recorded by the OPERA and MINOS collaborations HE]. When the 
Lagrangian includes vector self- interaction in the form of (fl2l) or (I29I) . the 
energy spectrum of a 1-dimensional neutrino beam (1281) includes additional 
term (13 7p . and the neutrino can move faster than light. Its energy spectrum 
is compatible with the power energy spectrum ([H]), and the neutrino velocity 
is associated with the observed data when the coupling constants satisfy 
constraints (fldj) and (H2l) . The Heisenberg model of self-interacting spin-1/2 
particle ( )32|) can also give satisfactory result under appropriate choice of the 
coupling constant (H2~j) 

Neutrino with scalar self- interaction (J33J) is always subluminal (16Tj) . 
Scalar-vector self-interaction can result in the group velocity f[8"2"j) above 
the speed of light but it is much smaller than the observed values dlj-©. 

Although the only model of vector self-interaction (|29|) is suitable for 
explanation of superluminal neutrino, it looks very promising because this 
effect is universal and does not depend on the medium or external fields. It 
should be emphasized that the neutrino mass is not important here, and the 
observed superluminal velocity can be produced by a massless particle. We 
have applied the mean-field approximation (138|) for obtaining the effective 
energy spectrum (122]) of a plane- wave solution (1351) . It is enough for estima- 
tion of the neutrino velocity. Exact analysis of equation (13T)j) may concern the 
problem of neutrino flavor oscillations that is the subject for further research. 

The author is grateful to Erwin Schmidt for discussions. 
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